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ABSTRACT 
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This  paper  presents  a  specialized  algorithm  for  the  transhipment 
along  a  single  line  problem.  The  problem  is  a  specially  structured 
network  flow  problem  for  which  the  basis  structure  is  such  that  a 
greedy  algorithm  can  be  employed  for  solution.  The  specialized 
algorithm  is  on  the  order  of  a  hundred  times  faster  than  the  primal 

^  ,  r  -  /  z~  /-  - 

simplex  method  on  a  graph.  r>  ' ■  "  ' 
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The  problem  of  transhipment  along  a  single  road  L23  may  be 
formulated  as  a  minimum  cost  network  -flow  problem. 


min  ex 

s.t.  Ax  =  r 


x  >  0 


where  A  is  a  node-arc  incidence  matrix  for  a  network  on  2n 
nodes  with  5n-4  arcs  as  given  in  Figure  1.  The  vectors  x,  c, 
and  r  are  respectively  the  vectors  of  decision  variables,  cost 
coefficients  and  requirements.  The  underlying  network  consists 
of  supply  points,  i,  (i  =  l,2,...,n)  each  with  supply 
supply  s^  and  n  demand  points,  j,  (j  -n+1,  n+2,  . ..,n+n)  each 
with  a  demand  d^  .  Arcs  of  the  form  (i,j) 

i  =*  l,2,...,n,  j  =  n+1 ,  n+2,  .  .  .  ,  n+n  are  referred  to  as 
transportati on  arcs  and  arcs  of  the  form  (i,j),  i  = 
n+1 , n+2, . . . , n+n ,  j  =  n+1 , n+2, . . . , n+n  are  referred  to  as 
transhipment  arcs.  Arcs  (i,n+i+l)  1 =1 , 2, . . . , n-1 ,  and  arcs 
< i , i + 1 )  i  =  n+1,  n+2, . . . , n+n- 1 ,  have  cost  f  (these  arcs  are 
called  forward  arcs)  ;  arcs  (i,n+i-l),  i  =  l,2,...,n-l  and  arcs 
(i,i  — 1),  i  =  n+1,  n+2,  ...,  n+n-1,  have  cost  (those  arcs  are 
called  back  arcs);  arcs  <i,n+i),  i=l,2,...,n  (called  direct 
arcs)  have  cost  d  .  There  is  an  added  stipulation  that 

d  f  +  b  (  1  ) 
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and  that  total  supply  be  equal  to  total  demand. 


The  problem  is  essentially  one  of  determining  how  the  flow 
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distributes  itself  over  the  arcs  of  the  network.  We  shall  see 
that  this  may  be  determined  by  asking  the  following  questions: 

(1)  When  are  transhipments  necessary?  This  essentially  asks 
when  does  there  have  to  be  flow  along  the  transhipment  arcs? 

(2)  When  is  shipment  over  the  forward  and  backward  arcs 
preferable  to  shipment  along  the  direct  arcs? 

In  answering  the  first  question  we  obtain  a  decomposition  of 
the  flow  into  flow  on  arcs  necessary  for  feasibility  and  flow 
satisfied  by  transportation  arcs  alone.  The  answer  to  the 
second  question  determines  the  distribution  of  the  rest  of  the 
flow  among  the  forward,  back  and  direct  arcs. 

In  order  to  answer  the  first  question,  we  examine  the  basis 
structure  for  the  problem!  The  number  of  basic  arcs  for  this 
problem  is  2n  -  1,  of  which  at  least  n  must  be  transportation 
arcs  and  at  most  n-1  can  be  transhipment  arcs.  Because  of  the 
conditions  (1)  and  (2),  transhipments  arcs  will  never  be 
admissible  for  basis  entry  in  primal  simplex  pivots.  As  such, 
the  only  flow  on  these  arcs  possible  is  flow  due  to  the 
suppl y/demand  structure.  That  is,  flow  necessary  for 
f easi bi 1 i ty. 

A  feasible  solution  which  pulls  back  flow  or  pushes  forward 
flow  along  a  transhipment  line  minimally  is  trivially  obtained. 
If  excess  supply  exists  at  some  point,  then  this  needs  to  be 
pushed  forward  along  the  transhipment  line.  If  excess  demand 
exists,  then  this  demand  must  be  met  by  pulling  flow  back:  along 
the  transhipment  line.  By  minimally,  is  meant  the  minimal 
amount  of  flow  required  to  be  along  transhipment,  arcs  to  ensure 
feasibility.  For  example,  consider  the  problem  in  Figure  2. 
This  problem  requires,  for  feasibility,  that  the  transhipment 
arcs  take  on  the  flow  given  in  Figure  3. 


Flow  in  arcs  to  ensure  feasibilit 
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After  the  determination  of  transhipment  flow,  the  supply  at 
node  1  is  equal  to  the  demand  at  node  n  +  i.  The  question  is  then 
whether  or  not  shipment  along  the  forward  and  backward  arcs  is 
cost  beneficial  when  compared  with  direct  shipment.  Note  that 
i  f 


d 

l 


f  + 

l 


b  .  , 

l 


then  direct  shipment  is  preferable. 


It  is  well  known  that  if  a  constant  is  added  to  the  costs  on 

all  arcs  which  emanate  from  a  node,  then  the  solution  to  the 

problem  remains  invariant.  Thus,  by  subtracting  d.  from  all 

arcs  which  emanate  from  node  i,  i  =  1 , 2, 3, . . . , n ,  we  have  a 

problem  in  which  all  direct  costs  are  equal  to  0  and  the  costs 

on  forward  and  back:  arcs  are  given  by  f  ’  and  b  7  where 

l  i 

f  7  =  f  -  d  i=l,2,..., n-1 

ill 

b  7  =  b .  -  d  _ ,  i=l,2,..., n - 1 

i  i  i+l 


For  two  pairs  of  suppl/  and  demand  nodes,  l,  n+i  and  i+l,n+i+l, 
the  question  of  whether  the  forward  and  backward  arcs  are  cost 
beneficial  is  answered  bv  e::amininq  the  quantity, z  =  t  .  7  +  b .  7  . 

l  i  l 


f  '  +  b  7 

l  l 


(.i  ,  then  arcs  <i,n+i  +  l)  and  <1  +  1,  n  +  i) 
may  have  flow  on  them. 


It  is  now  apparent  that  the  most  cost  effective  arcs  tor  the 
entire  problem  are  arcs  <k,n+k  +  l),  (k  +  l,n  +  k)  where  k  is  such 
t  hi  a  t 


l 
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If  all  the?  quantities,  z  are  >  0,  then,  direct  shipments  are 
optimal.  Thus,  by  ordering  the  quantities,  in  ascending 

order,  we  obtain  a  natural  greedy  method  f  or  the  solution  o-f 
tne  second  part  of  the  problem:  In  succession,  ship  as  much  as 
is  possible,  while  ensuring  -feasibility,  along  forward  and  back 
acrs  which  correspond  to  the  most  negative  z  .  When  there 
are  no  more  negative  z  left,  make  remaining  shipments  along 
the  direct  arcs. 

The  above  translates  into  the  following  algorithm  for  the 
solution  of  the  transhipment  along  a  line  problem:  The 

variables  back,  forward  and  direct  represent  the  flow  on  back 
forward  and  direct  arcs  r espect i vel y .  The  back  and  forward 
arcs  have  flow  which  represent  the  net  flow  being  shipped 
from  point  i+I  to  i  and  from  point  1  to  i  +  1,  respectively. 


Step_0  Initialize. 

(a)  forward (1 ) ,back !i ) .direct  (i )=0,  f 1 ag ( i )  =  off, 
i  =  1 , 2 ,  .  .  .  ,  n . 


Styp_i  Determine  flow  necessary  for  feasibility. 

Far  l  =  1  to  n-1 

(a)  Pull  back.  If  (supply(i)  <  demand (n+i ) ) 

back (i  )  -  demand (n+i ) -suppl y ( i ) 
supp  1  yu  +  1)  -demand  (n  +  i  +  1 )  -back  ( l  ) 
demand ( n+ i ) -demand ( n  +  i ) -bac k  ( l ) 

(b)  Push  forward.  If  .supply (i)  demand (rd  l  )  i 

f  orwar  d  i  )  '-supp  l  y  ( l  )  -demand  (n  +  i  ) 
demand  (n+  l  +1  )  -demand  m  +  i  +1 )  -forward  <  i  ) 
r-L',pp  i  y  (  i  )  -supp  i.<i  )  -+  or  ward  ( l  ) 

(c )  1  (  ( supp 1 v  ( l  i  -  demand (n  +  \  ) )  r  1  i p . 
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Step _ 2 _ Determi ne  rest  of  flow  on  -forward  and  back  arcs. 

(a)  Lot  I-  be  such  that 

c  ( I. )  -  mi  n  *>.  f  ( 1  )  -d  ( 1  )  +b  (  i  )  —  d  (  i  +1 )  ,  i -1 ,  .  .  ,  n-1 , 

•f  1  ag  <  l  )  -o-f  f  .  ./ 

( b  )  FI  ag  <  k  ) =on . 

KC)  1  f  1 2  <  k  >  <  0) 

Determine  maximum  -flow  which  can  be  put  on 
f orwara ( k ) , back  (k ) ,  and  make  this  allocation, 
i-f  (z  (k)  >  O)  go  to  step  3. 

Step_3  Determine  -flow  on  direct  arcs. 

For  all  i  such  that:  suppiyfi)  >  0,  di  rect  (i  )  -suppl  y  ( i  )  . 

It  is  straightforward  to  see  that  all  steps  of  the  algorithm 
are  1 i near  except  for  Step  2  which  requires  a  sort  and  a  search 
to  determine  the  flow.  Step  2c  involves  a  simple  computation 
of  the  minimum  of  the  supplies  at  point  k  and  point  k  +  1.  If 
there  is  a  series  of  points,  say  1 , 1 + 1 , . . . , 1 +t  such  that  r(l), 
c  (1+1 zU+t)  all  have  the  same  value,  then  there  may  be 
alternate  assignments  of  flow  for  these  points.  In  such  an 
event,  when  the  flow  for  s'.i-l)  or  c<l+t  +  l)  is  to  be  assiqned. 
Step  2c  involves  the  comparison  of  the  flows  on  the  forward  and 
back  arcs  of  points  1,1+1  and  1+2. 
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Implementation  of  the  Algorithm: 

.'he  algorithm  has  been  coded  in  standard  FORTRAN  and  been 
tested  on  randomly  generated  problems.  The  computational 
testing  is  summarized  in  Table  1.  For  purposes  of  comparison, 
the  problems  generated  were  also  solved  by  a  NETFLQ  C 1 J , a 
general  purpose  code  for  the  solution  of  minimum  cost  network 
flaw  problems.  The  test  problems  generated  had  supplies  and 
demands  drawn  from  a  uniform  distribution  (5,100)  and  costs 
distributed  uniformly  between  1  and  50.  The  computational 
times  are  very  stable  and  show  that  the  algorithmic 
specialisation  that  results  for  the  transhipment  along  a  line 
problem  is  on  the  order  of  a  hundred  times  faster  than  the 
unspec i ai i z ed  primal  simplex  algorithm  on  a  network. 


-J 
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Problem  Seed  Nodes 


1 1 1 1 1 

50 

12345 

50 

34367 

50 

33889 

50 

11229 

50 

15834 

1 00 

33890 

100 

67811 

100 

98123 

100 

58923 

1 00 

58923 

150 

93451 

150 

12345 

150 

1 1  1  1  1 

150 

78600 

150 

Optimal  Solution 


146514 

128559 

103761 

159346 

741488 

153026 

1105341 

626802 

618579 

1094640 

2285600 

1876195 

647033 

1344982 

996977 


Solution  Time* 


@Ceedy_Al g. 

NETTLO 

.  004 

.  129 

.  003 

.  190 

.  003 

.  122 

.  003 

.  136 

.  004 

.220 

.010 

.213 

.  010 

.685 

.010 

.559 

.009 

.  660 

.  01  1 

.571 

.020 

1 . 324 

.  020 

1.699 

.  020 

1 . 305 

.  019 

1.353 

.  019 

1.592 

♦CPU  seconds  on  the  Dual  Cyber  170/50  at  the  University  of 
Texas  at  Austin. 
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